Spin paramagnetism in d-wave superconductors 
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The Ginzburg-Landau equations are derived from the microscopic theory for clean layered su- 
perconductors with d x 2_ y 2 pairing symmetry, including the Pauli paramagnetism effect. The upper 
critical field H c i parallel to the c axis is calculated. A comparison with the experimental data for 
YBCO suggests that, relative to the orbital effect, the Pauli paramagnetism contribution to H C 2 
is significant. The reversible magnetization M in high magnetic fields is also calculated, showing 
strong temperature dependence of the slope dM/dH, as a consequence of the spin paramagnetism. 
A simple expression for the high temperature spin susceptibility is derived, in a good agreement 
with the Knight shift measurements on YBCO. 



PACS: 74.72.-h; 74.25.Ha 



I. INTRODUCTION 



^ • Among the most important properties of the cuprate superconductors are their magnetic properties, reflecting 
the highly anisotropic layered structure and unconventional pairing mechanism and symmetry.El Nowadays, the d- 
wave pairing symmetry in hole doped cuprates,0 and quasi-two-dimensional (2D) nature of superconductivity^ are 
well established. .Although the pairing mechanism is not known, the weak couplijig BCS model for 2D d-wave 
superconductivity,!!! and the corresponding Ginzburg-Landau (GL) approximation,U~B are remarkably successful. 
' The existence of superconducting phase in very high magnetic fields and layered structure of high-T c cuprates, makes 
^\ , the paramagnetic effect in the superconducting state much more important than in conventional superconductorsJj't£l 
The purpose of this paper is to derive GL equations for layered d-wave superconductors including the Pauli paramag- 
netism effect, giving simple analytical expressions for the upper critical field, magnetization, and spin susceptibility, 
suitable for comparison with experiments. 

The GLi equations for conventional (isotropic 3D s-wave) superconductors were first derived from BCS theory by 
Gor'kovO For 2D clean superconductors with d-wave pairing GL equations are derived by Ren et alp and extended by 
Won and Maki,[] and Shiraishi et al.fi to include the higher order derivative terms. However, in the above references 



the magnetic field influence on the electron spins has been neglected. In theji^jjlassical papers Maki and Tsuneto 



studied the effect of the Pauli paramagnetism in conventional superconductorsEj'Ej Their results mainly refer to the 
dirty limit. 

In many underdoped and overdoped cuprates, where the upper critical fields H C 2 and transition temperatures 
T c are relatively low, the upward, positive curvature in the temperature dependence has rbeen observed ,oL3 unlike 
the Werthamer-Helfand-Hohenberg theory prediction for the conventional superconductors Jl3 A number of theoretical 
models for this unusual behavior has,been proposed, involving, for example, the influence of the scattering by magnetic 
impurities and of inhomogeneities,II3 the presence of d xy pairing symmetry,ILa and the higher Landau levels effect .lla 

In the optimally doped cuprates, the temperature dependence of H C 2 parallel to the c axis appears to be qualitatively 
the same as in the conventional superconductors. However, high values of H C 2 and the irreversibility effects in 
optimally doped cuprates pose an obstacle to the study of the field induced transition to the normal statec3 The 
relevant magnetization measurements of H C 2 on high quality single crystals of YBapCiifiOy-a are performed by Welp 



et al. 
et 



21 



in a temperature interval of about 8 K below T c . Recently, Nakagawa et a/.E3c3 for field parallel and O'Brien 
2 3 for field perpendicular to the c axis, reported data from GHz transport measurements up to 150 T of H — T 
phase diagram for Y r Ba2Cu307_,5 thiiuiilms in the whole temperature range. 

Yang and Sondhijj and Won et alS3 studied theoretically the paramagnetic state of d x i_ y i superconductors, ne- 
glecting the coupling of the magnetic field to the orbital motion of electrons in the superconducting planes. For 



perpendicular field, measurements support their theory, strongly: suggesting that H C 2 is limited by the spin paramag- 
netism below a certain characteristic temperature T* ~ 0.85T c .u For parallel field the role of the spin paramagnetism 
should be clarified with regard to the orbital effect. ._. 

Magnetization measurementg-pn YBCO in intermediate fields of Sok et al.pB are in a good agreement with the 
GL-like Hao and Clem model.E£l In higher fields, where the magnetic phase diagram contains a vortex fluid,E2l one 
can presume that the diamagnetisnuaf the vortex fluid resembles closely to that of an ideal Abrikosov mixed state 
at temperatures not too close to T C .E3 In this regime, the measurements on cuprates should have the same slope in 
different fields, with strong temperature dependence due to the spin paramagnetism. _ . 

Electronic spin susceptibility measuremoats provide evidence about the pairing stateO Recent ESR Knight shift 
measurements on YBCO of Janossy et alSB clearly show the temperature dependence characteristic for 2D d-wave 
superconductivity, different from the BCS theory prediction for 3D s-wave superconductors.^] 

In Sec. II we derive the GL equations for clean 2D d x i_ y i superconductors in the magnetic field parallel to the c 
axis, extending the procedure used by Ren et al.u to include the effect of the spin paramagnetism. We calculate the 
upper critical field and discuss the importance of the spin paramagnetism regarding the experimental data for YBCO 
in Sec. III. Within the model of layered superconductors as a stack of identical 2D conducting planes, we calculate 
the reversible magnetization in the Abrikosov approximation, Sec. IV, and the spin susceptibility, Sec. V. Section VI 
contains the conclusion. 



II. GINZBURG-LANDAU EQUATIONS 

Including the interaction of the electron spins with the magnetic field, the four coupled Gor'kov equations are 

G n (x,x> n )- J dx"A u (x,x")^+(x",x', W „)=5(x-x'), (2.1) 



iw n - — (-iV + -A) 2 - li b H + fx 
2m c 



1 e 
iu) n - — (-iV + -A) 2 + n B H + /j, 
2m c 



G u (x,x>„)- ydx"A iT (x,x")^ T + i (x",x', W „)=5(x-x'), (2.2) 



1 e 
-ioj n - — (iV + - A) 2 + li b H + ii 
2m c 



f +(x,x',o; n ) + J dx"A+ (x,x")G n (x",x',a; n ) = 0, (2.3) 



1 e 

-iw n - — (iV + - A) 2 - n B H + n 
2m c 



F+(x,x', Wn ) + J dx"A+ (x,x")G u (x",x',u, n ) = 0. (2.4) 



Here, G aQ (x, x', w n ) and F^(x, x', uj n ) are normal and anomalous Green's functions respectively)^ u> n = (2n + \)-kT 
are the Matsubara frequencies (fcs = % = 1), fi is the chemical potential and [ib is the Bohr magneton. The magnetic 



field H = curlA is perpendicular to the superconducting plane. The order parameter A^, is defined by self-consistency 



equation 



A+ (x,x") = V(x - x")(*{(xi)*{(x"t)) = -V(x - x")T^F+(x,x",c 



(2.5) 



where V(x — x") is the effective pairing interaction. Af» (x, x"), Af|(x, x") and A||(x, x") are defined analogously. 



From Eqs. fl2.1| ) — (2.5), to the third order in A, we find 

A+(x,y) =F(x-y)T^{ f dx'dx"G TT (x',x,-^„)A+ (x',x") [G u (x",y,( 



dxidx 2 G u (x ,xi,w„)A iT (xi,x 2 ) j dx 3 dx4G TT (x 3 ,X2,-w„)A 1 t i (x3,X4)G u (x4,y,cj„) 
Here G CT is the Green function in the normal state 



G CT (x,x',u; n ) = ^2 



ik(x-x') 



iljJn - (JUbH - £k 



3 -if A(x)-(x-x') 



(2.6) 



(2.7) 



where a — ±1 corresponds to spin up and down. The classical approximation is used, taking account of the orbital 
effect by a change of phase, and e k = (k 2 /2m) — fi. 
In the notation of Ref. 5, we find 



A*(R,k)=5>(k-k')T]T 



f 



,2 _ q, ,2 i 



(ojI + e 2 .,) 3 4m 2 



4m (w 2 +e 2 ,) 2 



n 2 



A*(R,k')-^nk-k')T^ 



f 



K 2 +£ 2^|A*(R,kO| 2 A*(R,kO, 



(2.8) 



where Af^ = A^, = A*, R = (x + y)/2 is the center of mass coordinate, and II = —id/dR — 2eA/c. TheFourier 
transform is performed with respect to the relative coordinate r = x y, and (ttT) 2 (fisH) 2 is assumed.113 
For d x 2_ y 2 pairing, we take 



v(k-k>) = v a Ckl-kl)(K 2 -K 2 )> 



and 



k 2 

x "'yj 



A*(R,k) = A5(R)(fc 2 
where k is the unit vector in the direction of k. 



(2.9) 
(2.10) 



From Eq. (2J3), we obtain for T near T c the first GL equation 



AS(R) = iiV(0)VoA51n^ - ? f ] J T2 N (0)V 



1 ,2tt2a* , „2 u2a» i Aq| 2 Aq 



y F n 2 A* + visH'At 



(2.11) 



where tu c is the energy cut-off, C — 1.78 is the Euler constant, and 2D density of states for a given spin orientation 
iV(0) = m/2n. This equation is the same as in Ref. 5, with an additional term containing the-Bphr magneton /j,b- 
In the same approach, from the quantum mechanical equation for a single electron of spin sE3 



.11. [(V**)^ - **W] - —A*** - 2c/x B curl(4'*s*), 
2m mc 



(2.12) 



the superconducting current density is 

j s (x) = 7j^^2^2 j rfx l dx 2rfx 3 dx 4 A(xi,X2)A*(x3,X4)G'_ (J (x3,X2,-W„) 



(V x - V y ) [G -(x,xi,w n )G ff (x4,y,( 



cTfiB curl ^""^ 5 j (ixi(ix2(ix3(ix4A(xi, X2)A*(x3, X4) 

u n a=-l J 

x G a (x, xi , ui n ) G- a (x 3 , x 2 , - uj n ) G a (x 4 , x, w„ ) , 



(2.13) 



where a is the unit vector orthogonal to the superconducting plane, a = all/ H. The first term in Eq. ( |2.13| ) is the 
diamagnetic current jf, and the second term is due to the difference between the paramagnetic current associated 
with the electron spins in the normal and the superconducting state, j s / m . Note that the total current is j = 
j s + 2A r (0)/i 2 3 curlH, where the second term is the paramagnetic current in the normal state. Introducing the center 
of mass, performing the Fourier transform of A with respect to the relative coordinates, and after integrations similar 
as in Ref. 5, we get 



eT x ^ 

Js = r / 



dk 



A*(R, k)- 



mt 



iuj n - e k 
2k 



k' = -k 



(iuj n - £k) 
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(2.14) 



Here, in the spin param agnet ic effect is neglected, as well as the orbital effect in 

For d x 2_ y 2 pairing, Eq. ( 2.10 ), after integration-over k, for T near T c , and assuming (ttT) 2 3> (hbH) 2 , the second 
GL equation for the total supercurrent density isL^I 



7C(3)iV(0)mw 



327T 2 T2 
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rM 2 B curl(H|AS| 2 ^ 



(2.15) 



Finally, for a stack of identical 2D conducting planes in the magnetic field parallel to the c axis, the corresponding 
GL equations in the standard form are 



^ + [3\ip\ 2 ip- 



1 

2m 



IT 2 V> + vH 2 ip = 0, 



(2.16) 



j* = ~^>% - - ^A|^| 2 - 2 C r;curl(H|^| 2 ), (2.17) 



9R mc 

where 



167r 2 T 2 T-T c 
7((3)mv 2 F T c ' 



(2.18) 



and 



48tt 2 T 2 
7C(3)A^*(0)m 2 4 



(2.19) 



7? = — (2.20) 



Here iV*(0) = N(0)/6, S being the average spacing between the planes.cl 

In this model, the GL equations for a d-wave, clean and layered superconductor are of the same form as in the 
isotropic and clean 3D s-wave case. However, besides the factor 2/3 in parameters a, (3, 77, the quantities m,VF, N(0) 
refer to the effective mass, Fermi velocity and density of states in 2D conducting ab planes. Therefore, the coherence 
length is £ = £ a & = (l/2m|a|) 1 / 2 , the penetration depth A = \ a b = (mc 2 j3 /Iftire 2 ^) 1 / 2 , and the GL parameter 

V6cT c 

(2.21) 



ev 2 F V 7C(3)iV*(G) 



The free energy density corresponding to Eqs. ( 2.16 ) and ( 2.17 ) is 



F = F n + c#| 2 + ||V| 4 + ^|n>| 2 + + VH 2 M 2 . (2.22) 



III. UPPER CRITICAL FIELD 



Near the second order phase transition to the normal phase, from the linearized Eq. (2.16) 

aib + — U* 2 iP + rjH 2 ^ = 0, 
2m 

we obtain for the upper critical magnetic field parallel to the c axis 



(3.1) 



H C 2 = 



8c/i 



/ 256^ 2 cV B r 2 T 
7C(3)e 2 4 {1 T c } 



(3.2) 



The slope at T c is 



dH c2 
dT 



IGfcT^ (33) 



and the GL expression without the Pauli paramagnetism correction is simply 



H° c2 =T c \dH c2 /dT\ Tc (l- J). (3.4) 



We illustrate our results using the experimental data foe YBCO-tlrtJ Samples are in the clean limit, and for 
fields above IT the spin-orbit scattering can be neglected.B Taking the slope dH^JdT\T C — —1-9 T/K from the 
magnetization measurements on optimally doped YBCO (T c — 92 K) of Welp et al.f3 we find that at T — 0.7 T c the 
spin pair breaking lowers the critical field for 10%, Fi g. 1 (a). However, it is evident that a larger slope corresponds 
to the data of Nakagawa et alE2 We obtain with Eq. ( |3.2| ) a good fit of the experimental data for T > 0.5 T c , taking 
the slope dH c2 / dT\x c — —2.6 T/K and T c = 84.3 K, Fig. 1(b). For this slope the paramagnetic correction is —15% 
at T = 0.7 T c . 

At low temperatures, where GL theory is not applicable, H C 2(T) becomes saturated. From the data of Nakagawa 
et aZ.,B H c2 (0) m 110 T. The Won and Maki expression H c2 (0) = -0.63 T c dH c2 /dT\ Tc P relevant for 2D d-wave 
superconductors, for T c — 84.3 K and the slope —2.6 T/K gives H c2 (0) = 138 T. Therefore, the spin paramagnetism, 
not taken into account in Ref. 7, but relevant in this case because the Clogston limit is 155-I.^hould lower H c2 (0) 
for 20%. However, note that H c2 measurements strongly depend on experimental techniques .ElEll 

The Fermi velocities corresponding to dH c2 /dT\ Tc = -1.9 T/K, T c = 92 K and dH c2 /dT\ Tc = -2.6 T/K, T c = 84.3 
K are vf — 10.1 • 10 6 cm/s and Up = 8.3 • 10 6 cm/s, respectively. The latter value is closer to that obtained from 
independently measured Fermi energy and the effective mass,0 giving vf — (7.6 ± 0.9) • 10 6 cm/s. 

IV. REVERSIBLE MAGNETIZATION 



The GL Eqs. ( p. 16 ), (2.17) in the dimensionless form are 

.V 



and 



(*--a) 2 / = /-|/| 2 /-7^ 2 /, (4-1) 

K 



curlcurla - | jfa - — (/*V/ - /V/*) + 7 curl(h|/| 2 ), (4.2) 
2k 



where 

,2 



327riV*(0)/4 T 



3 -^-Y} ^ 



The d ime nsionless quantities are introduced with tJj = ^J—a/[3f, R = AR_2e^A/c = a, and 2e£AH/c = h.Ei Eqs 



(H 



(4.3) are of the same form as for isotropic 3D s-wave superconducto with different k, and 7, Eqs. (|2.2l|) and 



Following the Abrikosov approachJ3 in the vicinity of the upper critical field parallel to the c axis of the layered 
superconductor, we obtain 

If I 2 

h = he- l £j, (4.4) 

T7|2 = 1 h c2~he , . 

1/1 k '(1-1/2k' 2 )/3' 1 ' 

b = h e - (h c2 - h e )j^—^. (4.6) 

Here, h e is the external magnetic field, (3 — |/| 4 /|/| 2 , the induction b is the mean value of the microscopic magnetic 
field h, and 



The dimensionless expression for the upper critical field, Eq. (|3.2|), is 



y/l+4K 2 7 - 1 

27K 



(4.8) 



We emphasize that-in Eq. (4.6) magnetization due to the spin polarization of the normal electrons is not included. As 



in 3D s-wave case,E3 the effect of the Pauli paramagnetism on b is included simply by the scaling of K — ► k'. Therefore 
the vortex lattice is again equilateral triangular lattice with f3 = 1.16. The temperature dep end en ce o f k' leads to the 
characteristic variation of the magnetization, 47rm = b — h e , with temperature. From Eqs. ( |4.6| )-(4.8), for n 1, the 
slope in physical units 



dM 



1 + 4k 2 7 



dM 



dH e 47t/3(2k 2 - 1) 



dH e 



Tc 



1 



256n 2 c 2 ^ B T ( 
7C(3)e 2 4 



(1 



T 



(4.9) 



is linearly decreasing with temperature, instead of the temperature independent slope in the GL approach with- 
out paramagnetic correction. Due to higher T c and smaller vf, this effect is not negligible in high-T c cuprates, 
in contrast to the conventional superconductors. For dH C 2/dT\x c — —1.9 T/K, we obtain a strong variation 
{dM/dH e )/(dM/dH e )\ Tc = 1 + 1.4(1 - T/T c ), or 1 + 2.6(1 - T/T c ) for dH c2 /dT\ Tc = -2.6 T/K. 

In high-T c superconductors the influence of the vortex lattice melting, and the strong fluctuation effects in the 
vicinity of T c also affect dM/dH e temperature dependence.E3 However, because of th e str ong paramagnetic effect, 
the characteristic increase of dM/dH e with decreasing temperature, predicted by Eq. (4.9), should be dominant for 
H > 0.5 H C 2, and 0.5 < T/T c < 0.9. This effect may be compensated-pnly by the strong temperature dependence of 
the effective parameter k 2 within the extended GL theory approach. Ej 



V. SPIN SUSCEPTIBILITY 



The current density due to the spin paramagnetism is 



jf in (r) = ccurl (Mf(r) - M„(r)), 



(5.1) 



where Mjp m (r) and M„,(jj) = 2jV*(0)/ i^H arc the magnetization du e to t he spin polarization in the superconducting 
and in the normal state.tU From Eq. ( |5.l| ) and the last term in Eq. (2.17), 



m: 



8 M 2 



(5.2) 



In the weak magnetic field 



-a/(3, and the spin susceptibility \s = M s s v%n jH is 



;(1 



T 

T~r. 



(5.3) 



This simple GL expression for the spin susceptibility of rthe layered d-wave superconductor is in a very good agreement 
with the microscopic theory result of Won and Maki,Q in a large temperature range T > 0.6 T c . Analogously, the 
susceptibility in the GL approach for an isotropic 3D s-wave, as well as for 2D s-wave, superconductor is 



T 



(5.4) 



in the agreement with the Yosida result from BCS theory.E] Therefore, d-wave and s-wave symmetries of the order 
parameter lead to different slopes of the spin susceptibility curve in the high temperature regime. _ 

This is illustrated in Fig. 2, in comparison with the ESR Knight shift data for YBCO of Janossy et altB Small 
discrepancy of the experimental data and the theoretical predictions at high temperatures could be the consequence of 
the strong-coupling effectsJ^J or of the presence of a small s-wave component of the order parameter which is possible 
in YBCOLI 



VI. CONCLUDING REMARKS 



The GL equations for layered d-wave clean superconductor are derived including the Pauli paramagnetism effect. 
The parallel upper critical field and the reversible magnetization in the Abrikosov approximation are calculated. The 
results are of the same form as obtained by Maki and Tsuneto in isotropic 3D s-wave case, with changed k and 7. 
The expression for H C 2 gives an useful correction to the standard GL result, providing the correct determination of 
dH C 2/ '<IT\tc for high-T c superconductors from the fit of experimental data in the relatively large range of temperatures 
below T c . From the comparison with the experimental data of for YBa2Cu30j_5 thin films^j we obtain dH C 2/ '<1T\tc — 
—2.6 T/K, and conclude that the effect of Pauli paramagnetism on H c2 parallel to the c axis is significant in comparison 
to the orbital effect, —15% at T = 0.7 T c . At zero temperature, this effect should be greater, about —20%. The strong 
temperature dependence of the magnetization slope dM/dH e due to the spin paramagnetism influence, should be 
experimentally detectable in high field measurements. We have also derived the GL expression for the Knight shift, 
reflecting the pairing symmetry in accordance with the experimental data. 
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FIG. 1. 

Temperature d epen dence of the parallel upper critical field of YBCO. Solid curve: GL result with the paramagnetic 
correction, Eq. (3.2). The weak coupling microscopic theory result without the paramagnetic correction (dotted 
curve),Q and the corresponding GL result, Eq. (3.4), (dashftd-Jine) are shown for comparison. Solid circles represent 

ilJB (a), and Nakagawa et alSB (b). 
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FIG. 2. 

Te mperature dependence of the normalized spin susceptibility. Solid 
(5.4). Dotted curves: the weak coupling microscopic theory calculations, 
and bottom curves to 3D s-wave case. Solid circles: data of Janossy et ali 



GL approximation, Eqs. (5.3) and 
Top curves correspond to 2D d-wave, 
for YBCO. 
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